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ABSTRACT  £  K 


r  The  local  approximation  order  from  a  scale  (S^)  of  approximating 


functions  on  is  characterised  in  terms  of  the  linear  span  (and  its 

Fourier  transform)  of  the  finitely  many  compactly  supported  functions  *9  P 


whose  integer  translates  ,  j  t  Z®  ,  span  the  space  S  ■  S1  from 
which  the  scale  is  derived.  This  provides  a  correction  of  similar  results 
stated,  and  proved  in  part,  by  Strang  and  Fix. 
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SIGNIFICANCE  AND  EXPLANATION 

In  multivariate  approximation,  we  are  still  groping  our  way  toward  an 
understanding  of  just  what  makes  a  given  family  of  approximants  do  a  good  job.  As 
an  aid  to  further  understanding,  the  people  in  the  field  consider  simple  models  of 
suitable  approximants  and  try  to  establish  under  what  circumstances  such  model 
families  can  approximate  every  smooth  function  well.  One  such  model  is  that  of  a 
scale  (s^)  of  approximating  spaces  Sh  ,  with  each  obtained  from  some  fixed 

space  S  by  scaling,  i.e.,  by  a  contraction  of  the  underlying  independent  variable: 


f  s(»/h) 


f  s  y. 


,v._\  , 

'••.‘v  v 


This  idealizes  the  situation  of  a  finer  and  finer  mesh.  A  second  idealization  is 
that  the  space  S  be  the  same  everywhere,  i.e.,  that  S  be  invariant  under  integer 
shifts.  While  spaces  actually  used  in  practice  are  often  more  complicated,  it  is 
hoped  that  the  knowledge  arrived  at  under  these  sing>lifying  assumptions  will 
nevertheless  be  relevant,  i.e.,  give  a  feel  for  what  is  to  be  expected  in  practical 
situations. 

A  very  simple  model  of  an  effective  scale  (S^)  of  approximants  on  if  is  the 
following:  S  is  spanned  by  the  integer-translates  of  one  compactly  supported 
function  <0  ,  and  the  space  Sy,  is  obtained  from  it  by  scaling,  i.e., 

Sh  f  -  j)  >  . 

The  approximation  order  obtainable  from  such  a  scale  is,  by  definition,  the 
largest  k  so  that 

dist(g,  Sh)  -  0(hk)  for  all  smooth  g  , 
with  the  distance  measured  in  some  suitable  norm. 


»  • 


.•  Vv\v 


Around  1970,  several  authors  concerned  with  an  analysis  of  the  Finite  Element 
method  characterized  in  several  ways  the  functions  <0  whose  associated  scale 
provides  approximation  order  k  .  These  results  are  gathered  in  papers  by  Fix  and 
Strang.  The  latter  go  further,  though,  and  consider  also  the  case  of  a  scale  for 
which  S  is  the  space  of  translates  of  not  just  one,  but  of  several  conq>actly 
supported  functions,  a  case  of  more  direct  practical  interest  in  several  variables 
where  even  rather  simple  approximants  have  this  more  complex  structure.  Strang  and 
Fix  intend  to  characterize  controlled  approximation  order,  an  order  in  which  the 
size  of  the  coefficients  a(j)  of  the  approximant  is  to  be  constrained  by  the  size 
of  g  .  But  they  do  not  provide  a  complete  argument  for  their  characterization;  in 
fact,  their  characterization  is  wrong,  as  Jia  has  recently  shown. 

The  present  report  provides  several  characterizations  of  a  related  but  simpler 
concept,  that  of  local  approximation  order,  in  which  the  approximant  is  constrained 
to  have  nonzero  coefficients  only  for  <T( */h  -  j)  with  support  near  the  support 
of  g  .  This  provides  justification  for  recent  work  by  Dahmen  and  Micchelli  and  by 
Jia  concerning  the  approximation  order  of  spans  of  box  splines. 


it.'.  ’  t 
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The  responsibility  for  the  wordino  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 
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CONTROLLED  APPROXIMATION  AND  A  CHARACTERIZATION 
OP  THE  LOCAL  APPROXIMATION  ORDER 


C.  de  Boor  and  R ,-q.  Jia 

The  tarn  "control lad  approxlnatlon"  was  Introduced  in  1970  by  Strang,  in  [St] »  see 
also  [rs] .  It  concarna  approximations  of  tha  form 


T. 


iff* 


with  •  a  finita  collaction  of  functions  on  f  of  compact  support  and  ®*c  tha  function 
obtained  from  ®  by  convolution  with  soma  "sequence"  ci  R  — >  if1  ,  1  •  * •  > 

«**c  «*  7  tf>[»-j)c(j)  . 

If,  more  generally,  c  is  some  function  on  rf®  ,  we  will  still  just  write  <f>*c  instead  of 
the  correct  but  aura  complicated  «'*(c|am)  . 

The  function  u  to  be  approximated  Ilea  in  the  Sobolev  apace  h£(^>  with  norm 


,UVp*“  fj<k,tt,l.P 


where 


and 


lu 


I.  !■  T  •  •  lD**u1 

'3.P  1  M*3  P 


•u« 


lul 


0,p  ,tt,L jf) 


He  denote  by  w*  <|P)  the  subepace  of  Of  compactly  supported  functions. 

P.c  P 

The  approximations  are,  more  explicitly,  of  the  form 


with 


Tj,(x)  ««  T(x/h)/h*,p  . 

Concerning  the  degree  of  approximation  to  u  f  H*(^>  achievable  by  proper  choice  of  the 
wei^its  cj  ,  Strang  and  Pix  [SP»  Theorem  III  state  the  following  result.  In  ita  statement 
and  subsequent  analyeis,  the  normalised  multivariate  monomials  appear  often  enough  to 
deserve  an  abbreviation  of  their  own.  He  will  use 


to  stand  for  tha  function 

^  — >  S:  x  | — >  x°/al 

(and  will  use  standard  lulti-index  notation  throughout).  In  particular,  D®U®  ”  [1**®  » 

and  this  holds  even  when  6  $  a  since  then  £]  *"  ®  -  0  by  convention.  Further,  »j  will 

denote  the  collection  of  polynomials  on  ^  of  total  degree  <  J  •  Finally,  f  will 

*  »i£x 

denote  the  Fourier  transform  of  f  ,  i.e.,  f ( £)  t“  J  e  f(x)dx  ,  with  £x  the  scalar 

i r 

product. 

theorem  at.  Imt  •  be  a  finite  subset  of  «|“’(^)  '  Then  tha  following  are 
equivalent) 

(i)  there  exists  a  sequence  <ta)|a|<k  in  BP*n  ♦  which  satisfies 

(ia)  *Q(0)  -  1  ,  ♦„  -  0  on  2r*®\0  » 

(ib)  Bia_B  -  0  21  2«^  for  0  <  I  a|  <  k  . 

(ii)  There  exists  a  sequence  <♦(,)) a|<k  iS.  *P*n  *  whloh  satisfies 

H*  “  £e<«  Vo* 11 6  '  iss.  M  <*  • 

(iii)  There  exist  soma  finitely  supported  c,p  so  that  t  tpc^  satisfies 

f(0)  ^  0  ,  but  Daf  *  0  on  2s^\0  for  |a|  <  K  . 

(iv)  For  each  u  €  ,  there  exist  weights  cj*  so  that 

(iva)  lu  -  I  \*c^'8  2  *  const^  hk  'lul^  2  *  swO, ...,k-1i 

(ivb)  I  IcN,  <  const  lul  2  . 

<P  <p  2  2 

For  the  very  special  case  when  m  *  ♦♦  ■  1  ,  such  results  can  already  be  found  in 
[Sc] .  In  [SF] ,  the  special  case  when  ♦  consists  of  just  one  function  (but  ■  is 
arbitrary)  is  treated  first  (see  [SF;  Theorem  I])  and  completely.  However,  for  the  general 
case,  [SF)  only  give  a  proof  for  the  implication*  (i)  ~>  (Hi)  — >  (iv)  .  In  particular. 


the  validity  of  the  implication  (iv)  —  >  (i)  has  recently  been  questioned.  This  was 


finally  aattlad  by  [J]  who  shows  by  a  counterexample  that  (lv)  does  not  imply  (1)  in 
general. 

This  raises  the  question  of  how  to  modify  (iv)  to  obtain  something  equivalent  to 
(i).  This  ia  a  matter  of  changing  the  control  over  the  form  of  the  approximation  as 
expressed  by  (ivb) •  In  this  connection  it  is  very  useful  to  recall  that  [DM]  quote  Theorem 
ST  in  a  modified  form.  The  modification  of  importance  here  occurs  in  condition  (iv)  which 
they  require  to  hold  locally,  as  follows i 


(iv)'  For  each  u  €  there  exist  weights  so  that,  for  any  closed  domain 


G  C  , 


(iv.)'  '«*I,VclVGI  *  °°n*t  hktul)c,p(Brh<6)) 

i  wa  -  )>»  ‘  —•* 

holds  for  sosw  const  and  r  independent  of  h  ,  G  ,  u  . 

Here,  Bd(G)  i«  (x  €  if*  i  dist(x,G)  <  d}  . 

If  (iv)'  holds,  then,  with  a  reference  to  [St] ,  [DM)  say  that  •  provides  "controlled 
Lp-approxlmation  of  order  k  They  do  not  coaewnt  on  the  fact  that  (iv)'  ia  a 

strengthening  of  (iv),  and  refer  to  [SF]  for  a  proof  of  the  implication  (iv)'  “>  (i). 

As  it  turns  out,  (iv)*  does  indeed  imply  (i).  But  it  is  the  localness  rather  than  the 

control  that  does  the  job.  For  this  reason,  we  propose  here  to  abandon  the  notion  of 

"controlled  approximation  order"  in  favor  of  "local  approximation  order".  He  say  that  ♦ 
provides  "local  L^-approximation  of  order  k  "  in  case  the  following  condition  holds; 


(iv)"  For  each  u  €  W^d®)  there  exist  weights  ch.  so  that 
-  -  p  1,1  -  V  — — — 


(iva)" 


,u  ■  e<p  vci'P  < coMt  h’tiu,k,I 


c"(j)  “  0  whenever  dist(jh,  supp  u)  >  r 


holds  for 


const  and 


r  Independent  of  h  and  u  ,  with 


*  e  "  •  "  *  *  !  •  ‘m 

W*n 

e* 


■.  v»  *  • . 


. •  . 

.• 

*.*  '.***  .  *  . 

* ■  .'*  ^ 

V^V-V-'s 


%  V*  •.*  *  * 

•  -  -  »  •  .  '  .  * 

"■ 

.  VV' 

»  *  •  -  “  * 

*. V- .V  *- 


.**’.*»  *•  ■  • 


•  V-.-vy 


^  ,/yyy.vy. 


'J  ^  J  J  V'  .-  .•  .•  .-  V  V  -  ■'  •  .  -  -  -  •  .7 


<P^(x)  1“  If( x/h )  . 

It  ia  clear  that  (ivb)'  la  atrongar  than  both  (ivb)  and  (ivb)",  but  (ivb)  and  (ivb)"  are 
not  comparable.  Aa  we  are  about  to  ahow,  (iv)"  ia  the  right  aodification  of  (iv)  to  give 
equivalence  with  (1).  Thia  ahowa  that  alao  (iv)'  “”>  (i)  and  ao  validatea  the  veraion  of 
Theorem  ST  in  [DM] . 

Theorem  Let  *  be  a  finite  aubeet  of  Wp(C(*^>  •  Then  the  following  atateaenta  are 
equlvalenti 

( 1  • )  There  exlata  a  aeguence  (♦0|)|a|<jt  ill  *P*n  ♦  which  aatlaflea 

(fa)  *0(0)  -  1  ,  -  0  on  2«*b\0  t 

(fb)  E0  [-1D]%  -  0  on  2rt"\0  for  0  <  |a|  <  k  . 

p«i  — 

(2*)  There  exlata  a  aeguence  (+a)|a|<k  —  BP*n  •  euch  that 

ll“"  Va  W"8  €  ^  l“l  <*  * 

(3*)  There  exiat  acme  finitely  aupported  c ^  ao  that  t  »•  Eggf  i^c^  aatiaflea 

tl*  -  ♦*[]*  €  a|fll|.1  for  | <x|  <  k  . 

(4* )  for  all  p  €  [1,*»]  ,  •  provldea  local  L^-approxlmatlon  order  k  . 

(5*)  r°r  ,0—  P  €  (1«*1  >  ♦  provldea  local  ^.-approximation  order  k  . 


.  Condition  (1*)  differ*  from  (i)  in  that  the  latter  requires,  additionally. 


,b: 


that  E. .  [-ID]  4  .  *  0  at  0  .  Already  [DM)  prove  that  it  ia  poaaible  to  get  away  with 

P'a  ct-B 


the  weaker  condition  (1*). 

Condition  (2*)  ia,  offhand,  weakar  than  (ii),  but  an  inductive  argument  leada  from 
(2»)  to  (ii). 

Condition  O*)  aeeme  more  uaeful  to  ua  in  application*  than  the  (equivalent)  condition 
(ill),  w*  note  that  induction  give*  the  aeemingly  atronger  statement 

(3*)'  There  exiat  soaw  finitely  aupported  c^  ao  that  t  : ■  E^f  aatlaflea 

♦•[)*  -  []°  for  |a|  <  k  . 


4  - 


The  assarted  equivalence  between  (4*)  and  (5*)  shows  that  one  might  as  well  drop  the 
qualifier  "Lp-“  and  just  speak  of  the  local  approximation  order  provided  by  6  . 


.  ...# 


Proof  of  the  then  real.  While  the  main  point  of  this  note  is  the  implication  (5*)  “> 
(1»),  we  give  a  proof  of  the  entire  ia$lication  cycle  (1*)  ”>  (2»)  ««>  ...  «»>  ( 1  • ) .  The 
arguments  for  (1*)  “>  (2*)  “>  (3*)  are  adaptations  of  those  in  [SFl. 

( 1  • )  “>  (2*).  Observe  that  the  Fourier  transform  of  4(x  -  • )  U  8  is 
(ID] ®(e-i* .  Hence,  by  Poisson's  Summation  formula,  we  have 

♦*tl6(x)  -  UD18(e’U*  i(-5)))e.2l|j  -  V2’ijx(-iD]8‘ri(-2itj)  . 

tt  follows  that 

• ')  *.<.  W  l'tD1 
-  V”1’"1,..  tl'r 


the  last  equation  by  i1°). 


(2*)  “>  (3*).  With  y  i*  (k,...,k),  the  monomial  []Y  and  its  integer  shifts  span 
the  space  r^dl)*  ...  .  This  implies  that  there  exist  finitely  supported  sequences 


eg  so  that 


[]T*C-  -  n*-e  for  8  <  Y 


On  applying  DY_a  to  both  sides,  we  find  that  n°*Cg  “  l]°~®  ,  hence 

cB*[]a  ■  [l®*8  on  for  a,  B  <  Y  • 


*  !“  E|8|<k  VCB  * 

Then 

’  E|8|<kVc8*n<*  “  ^iBhkV11*’8  “  EB4«Vt]a’8  €  »“  *  *|«|-l  • 

the  inclusion  by  (2*). 

The  proof  of  (3*)  "“>  (4»)  follows  the  argument  for  [BH»  Corollary  to  Theorem  6] < 

This  leaves 


»  .  m  "  m  • 

•/»  /  e  * 


»...  .  Ji. 


•  • 


1  *  '  "  '  “  w  •  *  »■ 


.\v\vj  .vl' v\ 

'.'.VA"  .*  s'  •*. 
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(5*)  ">  (1*).  we  approximate  a  tensor  product  of  univariate  B-splines,  namely  the 


function 


with 


u(x)  !-  *  Mjt+^(x(v) ) 


1“  — j— ,  -~2~!  • «  •  ,  — J— J  (  *  “  t  J+  ,  t  f  R 


»W*> 

(see  [Sc]).  Since  u  f  for  any  p  f  f1,“)  ,  we  can  find  weights  ck  so  that  (iva)" 

and  (ivb)*  hold  (whatever  the  p  might  be ) •  Set 

“h  !’  T*>  <VC* 

and  consider  the  Fourier-Laplace  transform 

f(z)  t  e”ilxf(x)dx  ,  z  f  (f  , 

of  the  error  f  :■  u  -  Uw  •  Since  u  has  compact  support,  so  does  u^  by  (ivb)",  hence  so 
does  f  ,  uniformly  in  h  .  This  means  that  supp  f  lies  in  some  ball  Ba  of  finite 
radius  a  independently  of  h  .  Consequently 


a  p 


while,  by  (iva)". 


This  implies  that 


|f(z)|  *  | e  lzxf(x)dx(  <  e*,lBZ,const_lf l_ 
a 


If!  <  const  h  |u|  • 

P 


I f ( z) |  <  const  hx  for  *Ims<  <  const 


we  can  therefore  invoke  Cauchy's  formula  (see,  e.g. ,  [R] )  to  get  the  estimate 

(1)  ,D<lf  *.(«)  <  const  hk  . 

The  Fourier  transform  of  u  is  well-known  (see  (Sc])r  it  is 

viz)  -  Tf  CUtv))1**1 
v-1 

with  C(t)  :»  (sin  t/2)/(t/2)  Whittaker's  Cardinal  function.  From  this  we  deduce  that 

(2)  u(0)  -  1 
and 

(3)  lim  (O^u) (x/h)/hk  -  0  for  x  €  I?\0  , 

h*0 

A 

hence,  by  (1),  must  satisfy  corresponding  conditions. 

We  now  compute  u^  .  Since 


tm  .  J 


& 


_r. 


-6- 


we  find  that 


je‘izx(p(x/h  -  j )dx  -  h“  4>(h*)e"ihz:> 


uh(*)  -  E^  V(h*)vJ  „(*) 

with 

<ou)  '■ h" 1  j  c5<^*’ih8j  • 

Thua,  from  (1)  and  (2), 

(4)  lia  E  <p(0)  v*  (0)  -  1  . 

h*0  v 

Further 

[D]a(i<h*>v*y*))  -  Eg<a  J6l[D]6i<h«)  <-ih)  la'B 

witn 

v5.y(z)  8"  h°  £3  c5(3>  tiJYe"ih£j 

each  a  (2»/h) -periodic  function.  (Note  that,  for  Y  "  0  ,  this  agrees  with  the  earlier 


definition,  as  it  should. )  this  implies  that 

vJ  Y(2sj/h)  -  v*  y«»  , 

and  therefore  ( 1 )  and  ( 3 )  give 

lim  E^  Eg<a  Jal<-i)la-6llD]65><2»J>  v^a_g{0>/hk  -  0  , 
h+0 

hence,  for  |a|  <  k  , 

(5)  lim  E  £  l-iD]B5>(2wj)  v£  g(0)  -  0  for  j  t  ^\0  . 

h-»0  v  p  9 

A 

From  (4)  and  (S),  we  deduce  (1*)  as  follows.  By  (4),  we  cannot  have  <p(0)  -  0  for 


every  <P  t  ♦  .  Without  loss  of  generality,  we  can  therefore  assume  that,  for  some  K  S  ♦ 

A  A 

X(0)  -  1  while  *K0)  -  0  for  all  <f  t  4\x  .  In  particular,  this  implies  with  (4)  that 
1l'*h+0vx  0* 0 '  "  1  *  ltow  consider  the  space  S  of  all  vectors  w  “  (w^^ )  for  which 


lim  E 
h+0 


(0) 


0  . 


We  claim  that  si  contains  a  vector  w*  with  ■  1  .  Indeed,  if  w£  0  ■  0  for 


all  w'  t  Si  ,  then  (Si)i  -  8  would  contain  the  unit  vector  (6^6.^)  •  hence 
1i*h*0Vx  o(0)  “  0  would  follow . 


With  this,  define 


£  w'  to 
<P  Y 


Then 

vo)  -  vi.oi(0)  ■  w;,ox(o,  - 1  • 

Further,  from  (5), 

£<p  £e<a  wi,a-6[‘iDl^<2,,j>  “  °  for  >“•  <  k  and  1  € 

and  this  says  that 

Zg<al-iD]  -  0  on  2*s“\0  for  |  a|  <  k  , 

as  we  wanted  to  show.  | | | 


Remark*  In  contrast  to  [SF] ,  we  assume  4  only  to  lie  in  1^  .  For  this  reason,  we 
do  not  get  (iva),  i.e. ,  we  do  not  get  simultaneous  approximation  to  derivatives.  But  this 
is  easily  obtained  under  the  assumption  that  ♦  lie  in  an  appropriately  smoother  space, 
using  the  quasi -interpolant  constructed  in  the  proof  of  (3*)  —>  (4).  In  this  connection 
we  note  that,  using  (3*)  in  the  equivalent  formulation  (3*)',  such  a  quasi- interpolant 
for  u  takes  the  particularly  simple  form 


.F  • 


•  •. 
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